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The problem of non-interacting electrons on a square lattice subject to a random magnetic flux is 
mapped onto a one-dimensional model with infinitely many orbitals per site. Linking each orbital 
with 1) other orbitals maps the problem onto Wegner's TV-orbital model in the same limit, 

while the original problem corresponds to TV = 1. The exact solution for N = oo, the mean-field 
limit, is discussed and compared with numerical results. An outline of a l/7V-expansion is given. 



The so-called random flux model (RFM) of electrons 
on a square lattice with a random magnetic flux penetrat- 
ing each plaquette has been the subject of some debate 
lately. Of particular interest is the question whether or 
not there are extended states and/or localization transi- 
tions in such a system. Some investigations of this model 
have been motivated by the quantum Hall problem at half 
filling: After a transformation to composite fermions, the 
fermionic quasiparticles at half filling are subject only to 
a fiuctuating vector potential, but not to a mean mag- 
netic field. ^ In the presence of disorder, the random fiux 
model is believed to provide a good representation of this 
situation, due to a constraint in the composite fermion 
theory that ties the fiuctuation of the vector potential 
to those of the fermion number density.^ Similar models 
of fermions with a gauge interaction have been proposed 
in the context of high-T^ superconductivity.^ Finally, an 
electron moving in a random fiux also provides an inter- 
esting localization problem in its own right. 

No definite conclusion about the nature of the states 
at the band center was obtained from various numerical 
investigations of the RFM. They were based on studies of 
the participation ratio, the conductance distribution, the 
scaling of transmission properties, quantum diffusion and 
Chern numbers. The existence of extended states near 
the band center was suggested. '^^^ However, data that are 
consistent with complete localization,®'® or with a whole 
region of critical states"'^'^ were also obtained. Analyti- 
cally, Aronov et al.^^ have shown that perturbatively the 
electron dynamics is diffusive. From this, they drew the 
conclusion that the model is in the same universality class 
as a supersymmetric nonlinear sigma-model with unitary 
symmetry, which is known to have only localized states 
in two dimensions (2D). However, the work of Zhang and 
Arovas"'^^ suggests that this argument may be misleading, 
for a similar reason for which its analog in the case of a 



homogeneous magnetic field fails. These authors argue 
that a term describing the fiuctuations of a topological 
density must be added to the nonlinear sigma-model, 
and that the system exhibits a Kosterlitz-Thouless tran- 
sition from a region with power-law localization at the 
band center to one with exponential localization in the 
band tails. 

In this paper, we propose an alternative approach 
to the RFM. We show that by means of a gauge 
transformation^^ the model can be mapped onto a one- 
dimensional (ID) model with infinitely many orbitals per 
site that is reminiscent of Wegner's #-orbital model. "'^'^ 
A modification that links each orbital to 2N other or- 
bitals provides, in the mean-field limit N co, an exact 
mapping onto Wegner's model in the same limit, where 
a solution is known. ^'^ The original RFM corresponds to 
N = I. The conductivity in the mean-field limit is non- 
zero everywhere in the band. While this is certainly qual- 
itatively different from the behavior at N = 1, we expect 
the density of states (DOS) to be less strongly depen- 
dent on N. Indeed, our mean-field result for the DOS 
agrees reasonably well with numerical results for N = I. 
We also outline a l/#-expansion that will provide a new 
way to obtain information about the transport properties 
of the RFM. 

We start with the tight binding representation of the 
RFM, 

H = -tJ2e^'-.^\^)(y\ . (1) 
(x.y) 

Here the |x) are electron states on the sites x of a 
square lattice with unit lattice constant. The sum is 
over nearest neighbors only, and t > is the hopping 
matrix element. The phases O^^y are independent ran- 
dom variables apart from the hermiticity requirement 
6x,y = —6y,x- We assume that they are independently 
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distributed within the interval [— 7r,7r]. Next we recall 
that the phase disorder can be restricted to ID by means 
of a gauge transformation.^^ Let ej be the unit vector in 

j-direction (j = 1,2), and define ^^XjX+e = ^'x ''- Now 
consider the states |x) = \xi,X2) = \r,m). From now 
on we assume the system to have a length R in the r- 
direction, with periodic boundary conditions. In the re- 
direction we assume a finite width M. A local gauge 
transformation 

|r, m) ^ e"'''^'"* |r, m) , (2) 
with ipr,m = X]™=/ ^r^" transforms the Hamiltonian into 
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The (/)r,m are linear combinations of the 6i]m \ 

'i>m(r) = (t>r,m = m + <^r,m " fr + l.m- 



(3c) 



(4) 



Apart from hermiticity, they are therefore also indepen- 
dent random variables that obey the same distribution 
as the oi^rn^ . In particular. 



0. 



(5a) 



(5b) 



The gauge transformation has eliminated the disorder 
in the r-direction. In the absence of the disorder in the 
Hamiltonian, i.e. the random variables /™^'™ inEq. (3a), 
we would simply have a ID tight binding model that can 
be diagonalized by means of a Fourier transform from 
sites r to wavenumbers q. The corresponding Green's 
function is^® 



G(°)(g,z) = [z- 2tcosq]-\ 



(6) 



with z = uj + iO, the complex frequency or energy vari- 
able. With the random matrix elements /™^'™ , the model 
written in the form of Eqs. (3) can be interpreted as de- 
scribing a ID system with M orbitals per site r. We are 
interested in the behavior in the 2D limit R, M ^ co. 

The model as defined by Eqs. (3) is formally very sim- 
ilar to Wegner's site-diagonal #-orbital model"'^'* (with 



N = M). We have deliberately chosen our notation so 
as to make this obvious. A crucial difference, however, is 
that in the #-orbital model, Wegner assumed all of the 
with m ^ m' to be equally distributed, while in 
our case / is nonzero only if m and m' are nearest neigh- 
bors. Wegner's results therefore do not directly apply to 
the present case. 

In order to proceed, let us redefine the random vari- 
ables /, Eq. (3c), 



t 6r 



9n 



J't>„ 



,(r) 



(3c') 



where the 4'm,m'{i') are independent random phases with 
equal probability distributions for all m ^ m' that link 
points (r, m) and (r, m'), and the gm-m' determine the 
range of / in m-direction. The original random fiux 

- 6m, -l) 
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(1) 



M(6„ 



model corresponds to g^, 
in the limit M ^ co. 

We now define a new model by choosing = 1. Then 
the Hamiltonian is given by Eqs. (3), however with ran- 
dom variables 



f 



r.r' 



t 6r 



J't>„ 



,(r) 



(7a) 



that are independently distributed with a second mo- 
ment, 



I rm,m rn,n \ ,2 c S S S S 

\Jr r' J s s' I — ^r,r' ^r,s ^s,s' ^m,n' ^m' ,n- 



(7b) 



Equations (3) and (7) now desribe Wegner's unitary, site- 
diagonal model. This has been achieved by means of the 
mean-field construction of linking each orbital to every 
other orbital on the same site. 

For the model described by Eqs. (3) and (7) in the 
limit M CO, an exact solution is known for both the 
single-particle Green's function, and in particular for the 
density of states, and for the transport properties."'^'* This 
is due to the fact that, as a consequence of Eqs. (7), 
all terms in the perturbation expansion for the Green's 
function that contains crossed impurity lines vanish in 
the limit l/M 0. The remaining diagrams can be 
summed exactly. 

The resulting Green's function can be written. 



G(q, z) = [z — 2t cos q + S(z) 



(8) 



where S is the self energy. The exactly solvable limit 
mentioned above corresponds to the self-consistent Born 
approximation for the self energy becoming exact, i.e. S 
is given by,"''* 



E(z)=t2^G(°)(g,z + S(z)), 



(9) 



with G°(q, z) from Eq. (6). By performing the g-integral 
we obtain a quartic equation for the retarded self energy, 
^r{uj) = S(z = w + iO), 
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(10a) 



Here we have assumed t = I, i.e. we measure all ener- 
gies in units of half of the ID perfect crystal bandwidth. 
Equation (10a) is easily solved, and from its solution we 
obtain the DOS N(ui) as 



N{ui) = - ImSij(w). 

TT 



(10b) 



Figure 1 shows our mean-field result, Eqs. (10), for the 
DOS together with the result of numerical calculations 
for the model with = Qm^ (see the definition after 
Eq. (3c')), and for models with additional couplings be- 
tween orbitals (m, m ± 2) and a given site r of the ID 
chain corresponding to N = 2 (see Eq. (11) below). The 
latter were obtained by numerically diagonalizing sys- 
tems with the sizes 40^, 80^, 120^ (N = 1) and 40^ 
(N = 2). The accuracy of the results is better than 
3% {N = 1) roughly 10% (N = 2). 

The Lifshitz-bounds of the spectrum for N = I (cf. 
Eq. (1)) are ±4. From Eqs. (7) we see that the disor- 
der induced by the random phases is t, as large as 1/8 
of the unperturbed bandwidth. This strong disorder is 
refiected in the DOS: only a weak maximum in the vicin- 
ity of w = 2t is left over of the square root divergence 
in the unperturbed DOS at the perfect ID crystal band 
edge, and the band edge is shifted by more than 50%. 
Lifshitz tails are, as usual, suppressed in our mean-field 
model. However, as can be seen from the numerical data 
in Fig. 1, the tails are extremely small, even in the orginal 
model, and the overall qualitative features of the numer- 
ical result are well reproduced by the mean-field DOS. 

There are, however, also important differences between 
the mean-field result and the exact numerical DOS: In 
the numerical result for N = I, the total bandwidth is 
about 10% larger than in the mean-field result. This can 
be qualitatively understood by the following argument. 
For gm = gm\ and no phase disorder, the band edge 
would be at w = 4t, and the DOS would have a van Hove 
singularity in the center of the band. In the presence 
of phase disorder, the latter is smeared out due to the 
absence of translational invariance. Simultaneously, the 
width of the band is reduced, since the average kinetic 
energy corresponding to the m-direction decreases as a 
consequence of the randomization of the phases. In the 
mean-field model, the phases are effectively even more 
randomized so that the band width is further dimin- 
ished. This argument is supported by first numerical 
results for models with an increased number of couplings 
in the m-direction. Already when next-nearest neigh- 
bors are taken into account (Fig. 1), the width of the 
band is reduced to approximately the value predicted by 
the mean-field theory with N = co. (In the center of the 
band convergence with respect to the system size is not 
reached for the data shown in Fig. 1, were for N = 2 only 
systems with sizes up to 40^ have been considered). 



The net result of the strong phase disorder is a den- 
sity of states that bears no similarity any more with that 
of a two-dimensional ordered system, but rather resem- 
bles that of a one-dtmenstonal system with a very strong 
rounding of the van Hove singularities at the band edges. 
As can be seen from the figure, this qualitative feature is 
well captured by our mean-field theory. 

The conductivity in the mean-field limit is simply 
given by the square of the DOS.^'* This means that the 
conductivity is non-zero everywhere within the band, 
and hence all states are extended. This feature is cer- 
tainly characteristic of the mean-field character of the 
model. Information about the transport properties of 
the model with g^ = gm^ can therefore only be ontained 
by calculating corrections to the mean-field limit. 

In order to set up a systematic perturbation theory 
with the mean-field model as zeroth approximation, we 
define the function g^ in Eq. (3c') as 







if < |m| < # 
else 



(11) 



This choice couples 2N orbitals on each site. In the limit 
M ^ oo, N = M yields the mean-field model solved 
above, while N = 1 yields the original model as defined 
in Eqs. (3). It is now possible to set up a systematic 
expansion in powers of l/N , along the lines of Ref. 16. 
While this should yield only small, quantitative correc- 
tions to the DOS, we expect an expansion of the con- 
ductivity in powers of the small parameter l/N to show 
indications of localization, starting at 0(l/#^). Such a 
calculation will be rather different from the perturbation 
theory of Ref. 11, and especially the dependence of the 
results on the location within the band should be very 
interesting. The fact the the above numerical DOS with 
N = 2 agrees already quite well with the mean-field re- 
sult suggests also that very probably the expansion with 
respect to l/N will converge quite rapidly. The results 
of such an approach will be reported elsewhere. 

In summary, we have constructed a mean-field limit 
for the random fiux model. We have shown that the 
mean-field DOS reproduces the qualitative features of 
exact numerical results. A systematic method to expand 
about the mean-field limit has also been proposed. 
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FIG. 1. The density of states at TV = oo (solid line) com- 
pared to the numerical results for N = 1 (•) and N = 2 (o). 
The accuracy of the numerical data is 3% and about 10% for 
N = 1 and N = 2, respectively. 
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